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Abstract 



In this paper we consider a coupled hydrodynamical system which involves the 
Navier-Stokes equations for the velocity field and kinematic transport equations for the 
molecular orientation field. By applying the Chemin-Lerner's time-space estimates for 
the heat equation and the Fourier localization technique, we prove that when initial 
data belongs to the critical Besov spaces with negative-order, there exists a unique 
local solution, and this solution is global when initial data is small enough. As a 
corollary, we obtain existence of global self-similar solutions. In order to figure out the 
relation between the solution obtained here and weak solution of standard sense, we 
establish a stability result, which yields in a direct way that all global weak solutions 
associated with the same initial data must coincide with the solution obtained here, 
namely, weak-strong uniqueness holds. 

Keywords: Liquid crystal fiow; global existence; blow up; stability; weak-strong 
uniqueness 
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1 Introduction 

In this paper, we study the following hydrodynamical system modeling the flow of ne- 
matic liquid crystal in M": 



where u and P denote the velocity field and the pressure of the flow, respectively, d de- 
notes the (averaged) macroscopic/continuum molecule orientation field, J^, A,7 are positive 
constants, and /(d) is a Ginzburg-Landau approximation function. The notation Vd Vd 
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dtu - lyAu + u • Vu + VP = 
atd + u- Vd = 7(Ad-/(d)) 
div u = 0, 



A div (Vd0 Vd), 






(1.3) 
(1.4) 
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denotes the n x n matrix whose (i, j)-th entry is given by did ■ djd (1 < i,j < n). As in [5] 
and [14], we assume /(d) = for simphcity. Besides, since the size of the viscosity constants 
V, A and 7 do not play a special role in our discussion, we assume that they all equal to the 
unit. 

The above system (|l.ip - (|1.4p describes the time evolution of nematic liquid crystal 
materials (cf. |12j). Equation (jl.ip is the conservation of linear momentum (the force 
balance equation). Equation (|1.2p is the conservation of angular momentum, in which the 
left hand side represents the kinematic transport by the flow field, while the right hand 
side represents the internal relaxation due to the elastic energy. Finally, equation (jl.3p 
represents the incompressibility of the fluid. This system was first introduced by Lin [T2] 
as a simplified version of the liquid crystal model proposed by Ericksen in [3] and Leslie 
in [lOj . and retained most of the interesting mathematical properties of the liquid crystal 
model. Some basic results concerning the mathematical theory of this system were obtained 
by Lin and Liu in [14] and [H]. More precisely, in [14] they proved global existence of 
weak solutions by using the modified Galerkin method combined with some compactness 
argument. Moreover, they also proved global existence of strong solutions if the initial data 
is sufficiently small (or if the viscosity v is sufficiently large). In [15] they proved that 
the one-dimensional space-time Hausdorff measure of the singular set of "suitable" weak 
solutions is zero. Recently, by using the maximal regularity of Stokes equations and the 
parabolic equations, Hu and Wang [8^ proved global existence of strong solutions to the 
system (|l.ip - (|1.4p for initial data belonging to Besov spaces of positive-order under the 
smallness assumption. Here we also refer the reader to see [1], [6], [11], [13], [16], [17], [7], 
[To] and the references therein for more details of the physical background of this problem 
and some different models of similar equations. 

As in the work of Hu and Wang [S], we let = Vd. Then, taking the gradient of (II. 2p . 
noticing the facts that F CD F — F^F {F'^ denotes the transpose of F) and 

j=i J j=i J j=i J 

for all i, fc = 1, 2, • • • , n, and applying the Leray-Hopf projector P to eliminate the pressure 
P, we see that the system (|1.1I) - (|1.4I) can be reduced into the following system: 

5tu- Au = -Pu- Vu-Pdiv(F^F), (1.5) 
dtF - AF = -u-VF - FVu, (1.6) 
(u,F)|t=o = (uo,Fo), (1.7) 

where Fq = Vdo. Recall that P = / + V(-A)-Miv, i.e., P is the n x n matrix pseudo- 
differential operator in with the symbol {Sij — ^|j5-)"j^i, where I represents the unit 
operator and Sij is the Kronecker symbol. Later on we shall consider the Cauchy problem 
(fT31) -(fr71). 

The purpose of this paper is to prove global existence and stability of solutions to the 
problem ()1.5|) - (|1.7p in the critical Besov space Bp^l'^"'^^{W^) of negative-order. It is easy to 
verify that (11.51) and (|1.6p have the same scaling property as the Navier-Stokes equations 
(which are equations obtained by putting d = in (|l.ll) - (jl.3p '). namely, if (u, F) is a solution 
of (II. 5p and (|1.6p with initial data (uq, Fq), then for any 5 > 0, by letting 

us = Su{5x, 5^t), Fs = SF{6x, S'h), (1.8) 
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we see that (u^, F^) is a solution of (jl.Sp and (|1.6I) with initial data {5uo{Sx), 5Fo{Sx)). The 
so-called critical space for the equations (|1.5I) and (|1.6p is a function space of (u, F) such that 
the norm in it is invariant under the scaling (ll.Sp . The so-called self-similar solutions are 
solutions satisfying the scaling relation u{t,x) = us{t,x), F{t,x) = Fs{t,x) (for all 5 > 0, 
X e M" and t>0). Obviously, if (u,F) is a self-similar solution, then we must have 

Uo(A.t) = X-^Uo{x), FoiXx) - X-^Fo{x). 

Such initial data do not belong to any Lebesgue and Sobolev spaces due to their strong 
singularity at x = as well as slow decay as — > oo, however, they belong to some 
homogeneous Besov spaces with negative order. This is the reason why we study the problem 
p.5p - (jl.7|) in the critical Besov space of negative-order. By making use of Chemin-Lerner's 
time-space estimates of the heat equation and the Fourier localization technique, we shall 
prove that when initial data (uo,Jo) belongs to the critical Besov space Bp:J+"/P(R") for 
some suitable p and q, there exists a unique local solution, and this solution is global when 
initial data is small enough. As a corollary, we get existence of self-similar solutions; see 
Section 3. In order to figure out the relation between the solution obtained here and the 
weak solution studied by Lin and Liu in |14) . we shall prove a stability result, which yields 
in a direct way that all global weak solutions with the same initial data must coincide with 
solutions obtained here, which is called the weak-strong uniqueness. 

Our main results are as follows (for notations we refer the reader to see Section 2): 

Theorem 1.1. (Existence) Let n >2, 2 < p < 2n and 1 < r < oo. Suppose that (uo, Fq) G 
Sp.r^"^''(R") and div Uq = 0. Then there exists T > and a unique solution {u, F) of the 
Cauchy problem (|1.5p - (|1.7p such that 

{u,F)e n £«(o,r;B-^+"/p+2/9(R")). 

l<g<oo 

Moreover, there exists £ > such that if |j(uo, fo)||^-i+"/p ^ £; then the above assertion 

holds for T — oo, i.e., the solution {\i,F) is global. Furthermore, if (u, i^) and {u.,F) are 
two solutions of (|1.5p - (|1.7p with initial data (uo,-Fb) o.'^^d {\1q,Fq), respectively, div Uq = 
and div Uq — 0, then there exists a constant C > such that 

||(u - u,F - F)\\^^^^ ,^,^-i+^/p+2/q-^ < C\\{uq - Uq,Fo - ^o)||^-i+"/p- 

Remark 1.1. (i) If (uo,Fo) belongs to the closure of iS(M") in i3p^g^"^'°(]R"), we actually 
have {n,F) e C([0, T], B^^+"/p(R")). 

(ii) Since the appearance of the term FVu in (11.61) , we need the condition 2 <p < 2n. 
For the Navier-Stokes equations, similar results hold for all 2 < p < oo. 

Now by the standard procedure, based on the uniqueness of solutions, allows us to 
deduce the existence of self-similar solutions of (|1.5p - (ll.7p . 

Corollary 1.2. Let n > 2, n < p < 2n. Suppose that (uo,Fo) G Bp.l^"^^(R") and 
div Uo = 0, and furthermore, Uq and Fq are homogeneous functions with degree —1, i.e., 

uo{x) = Suo{Sx), Fo{x) = 6Fo{Sx). 

Then the global solution {u, F) constructed in Theorem \l.l\ is a self-similar solution. 
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Theorem 1.3. (Blow-up criterion) Under the hypotheses of Theorem ] 1.1[ we denote by T* 
the maximum existence time. IfT* < oo, then for any 2 < p < 2n and 2 < q < oo satisfying 
- + ->§, we have 

II("'^)IIl<!(0,T-;B-1,+ "''''+^/') = 

For initial data (uo,Fo) G Bp^q^"^'^ [W], by TheoremlLll there exists T > such that 
the system ^^-^J\ has a unique solution (u, F) G f\i<q<oo2,'i{0, T; Bp)-'^'''^^'^''' {W')). If 
(uojFq) is additionally in the space then it is not difficult to see that (u, F) is also 

a weak solution (in the standard sense, see [15 )■ A natural question is the following: Do all 
weak solutions coincide with the one we obtained in Theorem II . 11 '' In order to answer this 
question, we establish the following stability theorem. 

Theorem 1.4. (Stability) Let n > 2, 2 < p < oo , 2 < q < oo and ^ + | > 1. Assume 
that (uo, Fq) and (uq, Fq) be two vector fields in L^(M") such that div Uq = and div Uq = 

0, and (u,-F) and (u,_F) be two weak solutions associated with initial data (uo,-Fb) OLi^d 
(uo,i^o); respectively. Let w = u — ii, E = F — F. If we assume further that (u, _F) G 
L?(0, T; B;;J+"/''+^/''(M")), then for anyO<t<T, 

\\{w,E)\\l.+2 /*|l(Vw,V£;)||i.dr < ||(wo,i?o)|li. xexp (c T ||(u, F)|r ^^,.,,^,,,dr) , 

Jo ^ Jo ' 

where wq = Uq — Uq and Eg = Fq — Fq, and C > is a constant. 

It is clear that if Uq = Uq and Fq — Fq, then Theorem 11.41 implies that w = E = 0, 

1. e., u = u and F = F. Hence, we have the following weak-strong uniqueness result for the 
system ([T3)) - pr7|) . 

Corollary 1.5. (Weak-strong uniqueness) Let n > 2, 2 < p < oo, 2 < q < oo and 
^ -I- I > 1. Assume that (uo,Fo) G L^(R") (div Uq = 0), and {u,F) be a weak solution 
of the problem ()1.5p - (|1.7p with initial data (uo,i^o)- If we assume further that {u, F) G 
i^(0, T; Bp^g^"^''^^^^(M")), then all weak solutions associated with initial data (uo,-Fb) must 
coincide with {u,F) on the time interval [0,T). 

Remark 1.2. For initial data (uo,i^o) e ^p;^+"^^(M"), when 2<p<2ri, l<r<oo such 
that - + - > 1, by Theorem 1 1.1[ there exists 2 < o < oo such that - -|- 2 > i and the 

system ^^-^J) has a unique solution (u, F) G ^'(O, T; ^j;:J+"/^+^/''(R")). Hence, in this 
case, by CoroUarv 11.51 weak-strong uniqueness holds for the system p.5p - p.7p . 

Organization of the paper. In Section 2, we recall some basic facts about Littlewood-Paley 
decomposition and Besov spaces. In Section 3, we present the proof of Theorem 1 1.1[ which 
yields existence of global self-similar solutions. In Section 4, we prove Theorem 1 1.31 Section 
5 is devoted to the proof of Theorem 11.41 

2 Preliminaries 

We first recall some basic notions and preliminary results used in the proof of our main 
results. Let S{W) be the Schwartz space and 5'(R") be its dual. Given / G 5(M"), the 
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Fourier transform of it, ^(/) = /, is defined by 

Let Vi = {(, e M", 1^1 < |} and 2?2 = {C e M", | < |^| < §}. Choose two non-negative 
functions ^ 5(R") supported, respectively, in 2?i and 2?2 such that 

We denote (/'j(^) = (f>{2~^£^), h — J-^^cj) and h = T^^ip, where J-^^ is the inverse Fourier 
transform. Then the dyadic blocks Aj and Sj can be defined as follows: 

A,f = m'W)f ^ 2^"" / hi2^y)f{x - y)dy, 

5,/ = i^{2-W)f = V" f h{Vy)f{x - y)dy. 

JR" 

Here D = (£'i,£'2, • • • and = i^^S^;, = -!)• The set {Aj,5j}jez is called the 

Littlewood-Paley decomposition. Formally, = Sj — Sj^i is a frequency projection to the 
annulus {\^\ ^ 2^}, and Sj = J2k<j-i is a frequency projection to the ball {|^| < 2^}. 
For more details, please reader to P and [9]. Let Z(M") = {/ e 5(M") : ^"/(O) 0, Va G 
(N U {0})"}, and denote by Z'{W) the dual of it. 

Definition 2.1. Let s e M, {p,r) G [l,oo] x [l,oo], the homogeneous Besov space Bp^{M.") 
is defined by 

i?;,,(R") = {/eZ'(R"): 11/11^.^^ < oo}, 

where 

IIJII^^, = |(E,ez2^'^''ll^-/lllO'^'' ^"'^ l<^<oo> 
[supj-ez2^1|Aj/||LP forr^oo. 

Remark 2.1. The above definition does not depend on the choice of the couple (</>, '0). Recall 
that if either s < - or s = - and q~l^ then {Bt, „(R"), || • || ) is a Banach space. 

Let us now state some basic properties for the homogeneous Besov spaces. 

Lemma 2.2. (Bernstein's inequality !T!) Let k G Z+. There exists a constant C independent 
of f and j such that for all I < p < q < oo, the following estimate holds: 

supp/c m < 2^} ^ sup < C2^-^-+^"(i/^'-i/«)||/|U.. (2.1) 

|Q|=fc 

Lemma 2.3. ([T]) Let k G Z+ and \a\ = k for multi-index a. There exists a constant Ck 
depending only on k such that for all s ^M. and 1 < p, r < c», the following estimate holds: 

C^'WJWb^^^ < < CkWdVU.^^. (2.2) 
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We now recall the definition of the Chemin-Lerner space £''(0, T; ^(M")): 

Definition 2.4. Let s £ R, 1 < p, q,r < oo, and let T > be a fixed number, the space 
£9(0,r;B^ ,,(M")) IS defined by 

£«(0,r;i3;_,(M")) := {/ e 5'((0, T), Z'(M")) : ||/IU,(o,t;S=,,(m.)) < 

where 

ll/ll£9(o,T;B|^) = (X/^''"' II^J''^lli'(0,T:LP)) 



I/IIl, 



L'!(0,T;i3= 



Remarks 2.2. (i) We define the usual space L'(0, T; Bp ^(-1^")) equipped with the norm 

q/r s 1/q 

(ii) By Minkowski's inequality, it is readily to verify that 

ll/ll£<^(o,T;B= J < II/IIl<j(o.T;B^„) if 9<^' (2-3) 



lL5i(0,T 



J < ll/ll£,(o,T;i3|.j if ^ < 9- (2-4) 



The following product between functions will enable us to estimate nonlinear terms 
appeared in (|1.5p and (|1.6I) . 

Lemma 2.5. ([2], [H]) Let 1 <p, q, r, qi, q2 < oo, si, S2 < ^, Si + S2 > and ^ = ^ + ^- 

Then there exists a positive constant C depending only on si, S2,p, q, r, qi, q2 and n such that 

ll/sll^^^p ^.^=i+=2-f J - '-^II-^II£''i(o,T;b;;i,)II5||£<!2(o,T;s;2^)- (2-5) 

Notations: The product of Banach spaces X x y will be equipped with the usual norm 
\\if,9)\\xxy = \\f\\x+\\g\\y, and if A" = J, we use ||(/, <7)lk to denote by \\{f,g)\\x>cx. For 
two nxn matrixes A = {aij)f.j^i and B = we denote A : B — J2i j=i'^ijbij- 

Throughout the paper, C stands for a generic constant, and its value may change from line 
to line. 



3 Local and global existence of solution 

In this section we prove Theorem ll.il Let p and r be as in Theorem 1.1, i.e., 2 < p < 2n 
and 1 < r < oo, and let 1 < g < oo. We choose a number 2 < qi < 2q such that — + - > |. 

qi p Z 

For a constant T > to be specified later, we denote Xt = ii-?! (0, T; 5^;^^+"^''+^/''' (M")). 
In order to establish the desired estimates in Xt, let us recall the solvability of the Cauchy 
problem of the heat equation: 

(dtu- Au = f{x,t), a;eM", t>0, 
I u{x,0) = uo{x), X e W\ 
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Proposition 3.1. {^) Let s E M. and 1 < p,qi,r < oo, and let T > be a real number. 
Assume that uq G Bp ^{M") 
(|3.ip has a unique solution 



Assume that uq G B^,^(M") and f G £«i (0, T; 5^+^^''' ^(M")). Then the Cauchy problem 



ue n £''(o,r;B^+2/'?(K")). 



Moreover, there exists a constant C > depending only on n such that for any qi < q < oo, 

ll"ll£.(0,T;i3;+^/') ^ ^(ll^olls.^ + II/IU,i(o,T;B;+^/«i- = ))- (3.2) 

Besides, if uq belongs to the closure of S{W^) m B^ ,,(R"), then u e C([0, T), ,,(R")). 

We also recall an existence and uniqueness result for an abstract operator equation in a 
generic Banach space. For the proof we refer the reader to see Lemarie-Rieusset [5]. 

Proposition 3.2. (0) Let X be a Banach space and ^ : X x X ^ X is a bilinear bounded 
operator, \\ ■ \\x being the X-norm. Assume that for any ui,U2 G X, we have 



|B(ui,W2)|U < Co\\ui\\x\\u2\ 



Then for any y E X such that \\y\\x < e < ihe equation u = y + B(u, u) has a solution u 
in X . Moreover, this solution is the only one such that \\u\\x < 2e, and depends continuously 
on y in the following sense: if \\y\\x l£ u = y + B(m, m) and \\u\\x < 2e, then 

\\u-u\\x < 1 „4g(7j l^~^ll-^- 

Now for given (u,F) G Xt, we define by Q{u,F) ~ (u,F), where {u.,F) is a solution of 
the following linear equations: 

atu- Au = -Pu- Vu-Pdiv(i^'^F), (3.3) 
dtp - AF = -u • VF - FVu, (3.4) 
(u,F)|t=o - (uo,Fo). (3.5) 

Proposition 3.3. Let (u,F) G Xt- Then we have {u, F) G Xt- In addition, the following 
estimates hold: 

||u|U, < lle^'^UolU, + C(||u||^^ + WFWl^), (3.6) 
\\F\\xr < \\e'^Fo\\xr+C\\u\\xjF\\x^. (3.7) 

Proof. We prove only the results for u, it can be done analogous for F. By the Duhamel 
principle, p.3p can be transformed into the following equivalent integral equations: 



VL{t)=e'^UQ- j eS'-^^^P{u-Vu){T)dT- f e^'-^^^Pdiv{F^ F){T)dT. 
Jo Jo 

Since we have assumed 2 < p < 2n, 2 < < cxd and - + — > |, we can apply Lemma 



by choosing si = ^-^id S2 = — 2+^ + ^ and Lemma [2.31 to obtain that 
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Similarly, by Lemmas 12.31 and 12. 5[ 



||div(F ^)ll£,i/2(o,T;Bp,J+"^''+''^''l) - "^ll^ -^ll£<'l/2(0,T;Sp^+"'''+''^'l) 

Hence, from Proposition 13. II we know u g Xt- Moreover, by using the boundedness of P in 
the homogeneous Besov spaces, p. 81) and p.9p . we get 

<||e*^Uo|U.+C(||uf^^ + ||Ff^^). (3.10) 
This proves Proposition 13.31 □ 

The Proposition 13.31 implies that G is well-defined and maps Xt into itself. Moreover, 
from p.6p and (13.71) . we know there exists a constant Cq > such that for all (u, F) E Xt 
and (u,F) = Q{u,F), we have the following estimate: 

||(u,F)|U, < \\{e'''no,e'''F^)\U^+C4{n,F)\\l^. (3.11) 

Case 1. (The small initial data). Taking T — oo and denoting X = Xoa- From Proposition 
13.11 there exists a constant Ci such that we can rewrite p. lip as follows: 

||(u,^)|U<Ci||(uo,Fo)||^-iw.+Co||(u,i^)||2.. (3.12) 

Now if we choose £ > sufficiently small such that Ci||(uo, i^o)||^-i+"/p < e < that 
is to say, ||(uo, i^o)||^-i+"/p < ^ < ic^' ^'^^^ Proposition 13. 2[ the system (|1.5p ~ (|1.7p 
has a global solution. 

Case 2. (The large initial data). In this case we shall use the Fourier localization 
technique to obtain existence of local solution of the problem ()1.5p - ()1.7p . To this end, we 
split uo = uoi + U02 such that Uo(0 = uol{|j|>2"} + uol{|5|<2"} •= "01 + U02, where Ip 
represents the characteristic function on the domain V. Similarly, we split _Fo — Fqi + Fq2- 
Since 2 < p < 2n, by using the properties of the Besov spaces, there exists N G Z+ such 
that Ci||(uoi,i^oi)tl^-i+"/p < see that 

||(e*^Uo,e*^Fo)|U, < + ||(e*^uo2, e*^Fo2)|U,. (3.13) 

Applying the Bernstein's inequality, there exists a constant C2 such that 

||(e*'^Uo2, e*'^i^02)||A'T — ll(e*'^uo2, e*'^Fo2)||^,j^g 

<C2(2^)/<^MI(e*^uo2,e*^i^02)IU,,(o,^^^-i+./P) 
<C22(2^)/«iri/9H|(uo,^^o)|U-iwP. 

Hence, if we choose T small enough so that C22'-^^^/i^T^^'^^\\{uo,Fo)\\^-i+,^/p < \e, i.e., 

^^ic22l + (2A^)M||(uo,i^o)|U-wJ ' ^^'^^^ 
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then we have ||(e*^Uo2, e*^Fo2)||A'T ^ §■ This result together with p.l3p yield the fact that 
for such T defined by p.l4p . we have [Ke'^^Uo, e*'^i^o) IUt < £• By applying Proposition 13.21 
again, there exists a local solution to the system (jl.Sp (|1.7p . 

If (u,F) G Xt is a solution of the system (|1.5p - (|1.7p . then one can proceed the same 
way as the proof of Proposition 13.31 to obtain that 

Pu-Vu, Pdiv(F^i^), uVF, FVu e £«i/2(0,T;B-3+"/P+4/9i(]R")). 
Hence, for any ^ < q < oo, we have 

iu,F) e £9(0,T;Bp;i+"/''+2/9(M")). 

Moreover, if (uq, Fq) belongs to the closure of iS(M") in i3p,r^"^^(M"), then we have (u, F) e 

Finally, we consider the uniqueness of solution. Note that in Proposition l3.2l we obtained 
only a partial answer to the uniqueness problem of solution, i.e., in the closed ball B2e, the 
solution of (|1.5p - (|1.7p is unique. Now we intend to get rid of this restrictive condition. 

Let (u, F) and (u, F) be two solutions of (ll.5l) - (|1.7p in associated with initial data 
(uo,Jo) and (uo,Fo)j respectively. Set w = u — u and E ^ F — F. Then {w,E) satisfies 
the following equations: 

- Aw + w -Vu + u-Vw + \7 ■ {E^F) + V • {F^ E) = 0, 
dtE-AE + wVF + VL-VE + EVu + #Vw = 0, 
w{x,0) = Wo (a;) = Uo(a;) - Uo(x), E{x,0) = Eoix) = Fo{x) - Fo{x). 

As the proof of Proposition 13. 3[ we can prove that 

llw • Vu + u ■ Vw + V • (E'^F) + V • {F^E)\\x^ 

< c{\\u\u + + \\F\U + ||i^lU.)||(w,£;)iu, 

and 

||w • Vi^ + u • V£; + EVu + F\/w\\xt 

< C(||u|U, + llulU. + IIFIU, + l|/'|U,)l|(w,i?)||;,,. 
Hence, by Proposition 13.11 we get 

ll(w,£;)|U, <Ci||(wo,£;o)|ls-i+"/P 

+ Co{\\n\U^ + \\u\Ur + \\F\U + \\F\U^)\\{w,E)\U,. 

Denoting M{T) := Co{\\u\\xt + \\u\\xt + \\F\\xt + II^IUt)- By the Lebesgue dominated 
convergence theorem, we know that M{T) is a continuous nondecreasing function vanishing 
at zero. Hence, if we choose Ti sufficiently small such that M{Ti) < then 

||(w,£;)|U, <2Ci||(wo,-Bo)IU-i+"/P. (3.15) 

Repeating the above procedure to the interval [0,Ti), [Ti,2Ti), . . . enables us to conclude 
that there exists a constant C such that 

ll(w,£;)|U^ <c||(wo,-e;o)|Ib-i+"/p. (s.ie) 

This implies the uniqueness result immediately. 
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4 The proof of Theorem 11.31 



Let 2<p<2n, l<r<oo, 2<g<oo such that + 2 > 3 Assume that 

P Q ^ 

||(u, -F')||£,(Q j,.g-i+n./p+2/q^ < oo . By thc embedding relation (|2.4p and the proof of Theorem 
ll.ll we see that 

ll(u,-P)lli«>(0,T;i3-^+"/'') ^ ll("'-^)IU~(O.T;B-;+"/'') 

< Ci||(uo,Fo)||^-l+„/p +Co||(u,F)||2^^^^^.^_i + „/,+,/^^ =M<0O. 

It suffices to prove that if \\{u,F)\\^^^^ ^^-i+„/p+2/q^ < oo, then T* > T. In other words, 
if T* < oo, then \\{u, F)\\ ^^^^ = oo. To this end, for any t £ [0,T), we 

take {u{x,t), F{x,t)) as a new initial data of the problem (ll.5p - (|1.7l) . and split u{x,t) — 
Ui(x, t) + U2{x, t) such that 

u(^,0 = iil|4|>2iv(^,t) + ul|^|<2iv(^,t) := ur(^,t) +u5(^,i). 

Similarly, we split F{x,t) :— Fi{x,t) + F2{x,t). Since 2 < p < 2n, by using the properties 
of the Besov spaces, there exists a sufficiently large constant N E N such that 

Ci||(ui(i),Fi(i))|U-;+./. <|. (4.1) 
On the other hand, if we choose T > t such that 

then we can obtain ||(e*'^U2, e*'^F2) H^'t+r^ < §■ This result together with (|4.ip . by Propo- 
sition 13.21 yield that there exists a constant depending only on e and M such that for 
any t g [0,T), the problem (|1.5p - (|1.7p has a solution on the time interval [t,t + T^). By 
the uniqueness we know that all solutions obtained in this way are equal in their common 
existence interval, so that the solution can be extended to the time interval [0,T + Te). That 
is to say T* > T, we complete the proof of Theorem [131 



5 StabiUty and weak-strong uniqueness 

The aim of this section is to prove Theorem 11.41 Let us recall the definition of weak 
solutions to the system (|1.5p - (jl.7p . 

Definition 5.1. The vector-valued function (u, F) is called a weak solution of (|1.5|) - (|1.7|) 

on R" X (0, T) if it satisfies the following conditions: 

(1) (u,F) e L°°(0,T;L2(Kn)) nL2(o,T;iJi(R")) := {WS), where H^{W) = ^l aCK") is 
the usual homogeneous Sobolev space. 

(2) (u,_F) satisfies the system (|1.5p - (|1.7p in the distributional sense, i.e., divu — Q in the 
distributional sense and for all v G C^(M" x (0,T)) and G G C^(M" x (0,r)) with 
divv = 0, we have 

udt'vdxdt — I I Vu : Vvdxdt / u • Vu • vdxdt 



JE" JO JM" JO 

T 

F'^F : Vvdxdt 
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and 

i-T I- 

F : dtGdxdt - / I \7F : VGdxdt - / I u \/F : Gdxdt 
) " 

i-T 

FVu : Gdxdt. 

(3) The following energy inequality holds: 

{\u{t)f + \Fit)\^)dx + 2ff {\Vuf + \VF\^)dxdT<f {\no\' + \Fof)dx. 

i" Jo JR" Jr" 

Remark 5.1. Formally, taking v = u and G ~ F, and adding them together, we get 
~f {\n{t)\^ + \F{t)\^)dx+ f (iVup + iVFHdx = 0, 

which implies the above energy inequality. Here we have used the fact AB : C ^ A : CB^ = 
B : A^G for any three n x n matrixes A, B and C, 

Let (uo, -Fb) e L2(M"), (uo, Fq) G L'^{W), and we denote by (u, F) and (u, F) two weak 
solutions in the space WiS associated with the initial data (uq, Fq) and (uq, Fq), respectively. 
Assume that {u,F) G i«(0, T; B;;g+"^^+^^'^(R")), where 2<p<ooand2<(7<cxD 
satisfying ^ + | > 1. Obviously, the above energy inequality yields that 



\\m),Fmh+2 / ||(Vu(r),VF(T))l|i.dr < ||(uo, Fo)||i.; (5.1) 
Jo 

||(u(i),F(i))||i.+2 r||(Vu(r),VF(r))||i.dr < ||(uo,i^o)|li.. (5.2) 



Let w ^ u — u, E = F — F , wq = Uo — Uq and Eq — Fq ~ Fq. To prove Theorem 1 1.4[ it 
suffices to prove the following result: 



Proposition 5.2. Under the hypotheses of Theorem \1.4\ we have 

U^it),Em\h+2 r||(Vw(r),Vi?(r))||i.dr 
Jo 

< ||(wo,£;o)||i. xexp(c / ||(u(r),i^(T))f ^,^„/,+,,,dT). (5.3) 
Note that by jEI]) and (|5?2|) . 

||(w(i),i?(t))||i.+2 r||(Vw(r),V£;(r))||i.dr=||(u(t),i^(i))||i. + |l(ii(i),nO)lli^ 
Jo 

+ 2 / ||(Vu(T),VF(r))||i.dT + 2 / ||(Vii(r),VF(r))|ji.dT-2(u(i)|ii(t)) 
Jo Jo 

- 2(F(t)|F(i)) -4 / (Vu(r)|Vu(r))dT-4 / (VF(r)|VF(T))dT 
Jo Jo 

< \\{uo,Fo)\\l. + ||(uo,Fo)|li. -2(u(t)|u(t)) -2(F(0|F(i)) 

-if (Vu(T)|Vu(T))dT-4 /" (Vi^(T)|V/'(T))dr. 
Jo Jo 
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Here (•!•) denotes the scalar product in In order to prove Proposition 15.21 we need 

to introduce the foUowing lemma. 

Lemma 5.3. Under the hypothese of Theorem \1.4\ the following equality holds for all t <T , 

(u(t)|u(i)) + {F{t)\F{t))+2 f (Vu(T)|Vii(T))rfr + 2 f {W F {t)\\7 F {r}) dr 

Jo Jo 

(uo|uo) + (Fo|Fo) - /" [ yv ■ Vw ■ udxdT + [ [ F^ F : Vudxdr 

Jo Js." Jo ./r" 

+ [ [ F^F : VudxdT - [ [ Vu : (F^F + F^F)dxdT 
Jo Jr" Jo Jr" 

yv-VF : Fdxdr + [ [ F : FVMvdxdr. (5.4) 



We shall use Lemma 1.1 in [S^ to prove Lemma [ 

Lemma 5.4. ([5]) Let n>2,2<p<oo and 2 < q < oo such that ^ + | > 1- Then fo 
every T > 0, the trilinear form 

(u, V, w) G X W5 X L'?(0, T; Bp;i+"/P+2/«(M")) ^ ( f u • Vv wdxdt 

Jo Jr" 

is continuous. In particular, the following estimate holds: 
/ u • Vv • wdxdt 

<'^\ML'iO,T-L^)\\^'^\\lHo'^^^^^^ 
+ \\^^\\l^O.T-,L^)\M1^1(o,T;L^)\\^"^\\l^O^^^^ 

+ Ml'{o,T:L--)\\^'^\\Iho'^^^^^^ (5-5) 

Remark 5.2. Checking the proof of Lemma 1.1 in 5 in detail, we find that the special 
structure of the Navier-Stokes equations (div u = 0) was not used, and (|5.5p holds in both 
scalar and vector cases. 

The proof of Lemma 15.31 

Let us consider two smooth sequences of ({u„},{F„}) (div u„ ~ 0) and ({u„},{F„}) 
(div u„ = 0) such that 

flim™(u„,F„) = (u,F) in L^O,T; H\W')), 

\ lim„^oc(u„,i^«) = {n,F) weakly-star in L°°(0, T; ^^(R"))^ 



and 



'lim™(u„,F„) = (u,F) in ^^(o, T; iji(R")) n L^O, T; BpJ+"/^+'/^(M")), 



[lim„^o,(u„,F„) = (u,F) weakly-star in {0,T; L^{W')). 
We split the proof into the following two steps. 
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step 1. Taking the scalar product with u„ and u„ of the equation (|1.5p on u and u 

respectively, after integration in time and integration by parts in the space variables, we get 

^ ((a,u|u„) + (Vu|Vfi„) + (u • Vu|u„) + (div(F^F)|u„))dT = (5.6) 



and 



^ ((5,ii|u„) + (Vii|Vu„) + (u- Vii|u„) + (div(F^i^)|u„))dT = 0. (5.7) 

Note that we have assumed that both Vu„ and Vu„ converge in L^{0,T; L^(M")) towards 
Vu and Vu respectively, it is obvious that 

lim f / (Vu|Vu„)dr + / (Vu|Vu„)dT) = 2 / (Vu|Vu)dr. (5.8) 
Since u„ converges to u in L''(0, T; i?p^g^"^^^^^''(M")), by Lemma [5^ one obtains that 



lim / (u • Vu|u„)(ir = / (u • Vu|u)(ir. 

"^°°Jo Jo 



/o Jo 
Due to the fact div u = 0, this yields, by Lemma [5.41 again. 



(5.9) 



(5.10) 



lim / (u • Vu|u„)(ir = — lim / (u • Vu„|u)(iT = — / (u • Vu|u)dr. 

Applying (|5.5p . it is also obvious that the following two results hold: 

lim / (div(F^i^)|u„)dr = - lim / (F^F|Vu„)dr 

= - I {F^F\Vu)dT^ [ {div{F^ F)\u)dT. (5.11 
Jo Jo 



and 



lim / (div(F^F)|u„)dT = lim / ( {d^. F^ F + F^ d^^. F)\un)dT 

= / iy^idx^F^F + F^d,^F)\u)dT = / {div{F^ F)\u)dT. (5.12) 
Jo ^ Jo 

Since dtu = Au — Pu ■ Vu — Pdiv(F-^F) holds in the sense of distribution, the estimates 
(|5.8[) - (j5.12p imply in particular that 

lim / (9^u|u„)dT = - lim / f (Vu|Vu„) + (u • Vu|u„) + (div(F^#)|u„)')dr 

= - J ((Vu|Vu) + (u- Vu|u) + (div(F^F)|u))dT 

— / {driL\u)dT. 
Jo 
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Similarly, we have 



lim / {drU\Un)dT — / {dru\u)dT. 



Putting these estimates together, and noticing that 

/ (a,u|u) + (9,u|u)dr = (u(i)|u(t)) - (uoluo) 
Jo 



(5.13) 



and 



we find 



and 



^ ((u ■ Vu|u) - (u • Vu|u)) dr = ^ (w Vw|u)dr, (5.14) 

{u{t)\u{t))+2 I (Vu|Vu)(iT = (uo|uo) - j (w • Vw|u)dT 
Jo Jo 

+ f f F^F : VudxdT + f f F^F : Vudxdr. (5.15) 

Jo Jk" Jq Jr" 

Step 2. In this step we derive the estimate for F and F. Proceeding the same way as 
I and (|5.7p . we obtain that 

^ (^{drF\Fr,) + (VF|VK) + (u • VF|F„) + (FVu|F„))dr = (5.16) 

^ ((9,#|i^„) + (V#|VF„) + (u. V#|^„) + (^Vu|F„))dT = 0. (5.17) 



Since we have assumed that both VF„ and VF„ respectively converge to VF and VF in 



1,2(0, r;L2(R")) and F„ converge to in Li{0,T- B~1^''^p^^^\MJ')). Then by using ^ 
and the fact div u = 0, it is clear that 



(5.18) 



lim f / (Vi^|VF„)dT+ / (V#|VF„)dT) = 2 / {VF\VF)dT, 

lim / (u- V/'|i^„)dT = / {u-VF\F)dT, (5.19) 
Jo Jo 

lim / (u ■ VF|F:„)dr = - lim / {u ■ VFn\F)dT = - [ {u-VF\F)dT, (5.20) 

n-^cx, _/q n^oo J^ Jq 

lim / {FVu\Fn)dT = / (FVu|F)dT, (5.21) 
"^°°Jo Jo 

lim / {FVu\Fn)dT = / (FVu|F)dT. (5.22) 

"^°oJo Jo 

The last identity holds because VF„ converge to V/" in L2(0, T; L2(M")) and {/'n} is 
bounded in L°°{0,T; L'^{R")), which was ensured by the Banach-Steinhaus theorem due 



and 
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to Fn weakly-star converge to F in L°°{0,T; L'^{M.'^)). Hence, combining the estimates 
(|5.18l) " (|5.22p . as derivation of estimate (|5.13p and (I5.14p . we have 



lim / {{drF\Fn)dT= [ {drF\F)dT (5.23) 



and 



hm / {drF\Fn)dT= [ {drF\F)dT. (5.24) 



It is obvious that 



/ {drF\F) + {drF\F)dT = {F{t)\F{t)) - {Fo\Fo). 
Jo 

Moreover, since div u = 0, we have 

/ [ (u • V# : + u • VF : f) dxdr = /" [ u- V{F : F)dxdT = 
Jo Jr" ^ 'Jo Jr" 

and 

FVu : F + F : FVu = Vu : {F^ F + F^F). 
These two facts imply that 

j (u-VF : F + U- VF : F + FVu : F + FVu : F^dxdr 

= [ [ (vu : (F^F + F^F) + (u - u) • VF : F - F : FV(u - u)) dxdr. 
Jo Jr" ^ ' 

Finally, putting all above estimates together yield 

(F(i)|F(i)) +2 / (VF|VF)dT = (FolFo) 

~ / / ' ^^^^ ^ ^^^^ +w -"^F : F - F : FVw) dxdr. (5.25) 

It is clear that (|5.4p follows from (|5.15p and (|5.25p . This completes the proof of Lemma 
The proof of Proposition 15.21 
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Note that J* Jg^„ w • VF : Fdxdr = since div w = 0, then by Lemma we get 

\\{Mt),Em\h+2 f ||(Vw(r), V£;(r))||i.dr < ||(uo, i^o)|li. + ll(£io, i^o)|li. 

-2(uo|uo)-2(Fo|Fo) + 2 / / w Vw • uda:rfr - 2 / / F : Vudxdr 

Jq Jm" Jq -/r" 

2 I f F^F : VudxdT + 2 f I Vu : {F^F + F^F)dxdT 



w ■ VF : Fdxdr ~ 2 



F : FVwdxdr 



<\\{wo,Eo)\\l2 +2 / wVw-udxdT-2 / / E^E-.Wudxdr 

Jo Jm" Jo Jr" 

w VF : Fdxdr + 2 f f E'^F : Vwdxdr. 

/O JK" Jo JR" 

Using the similar argument as the proof of Lemma [5T4l (see [5]), we obtain that 



- 2 



w • Vw-udxdr 



-l + Tl/p + 2/, 



<C /*||wf//||Vwf^7^/^||u 
Jo 

<\ I ||Vw||i.dr + C / ||w||i.||ur^,^„,,+,/,dT; 
^ Jo Jo P ' 

- [ I div(F^F) • udxdr 

Jo JK" 



F'^F -.Wudxdr 



< C / ||F||^/«||VF||^;'/'||u|L-i+„/P+v,rfr 
Jo 

< I f \\VE\\l.dr + C TllFlli^lluf 
^ Jo Jo P-'' 

n(w «) F) • V Fdxdr 

<l r||(Vw,VF)||i.dr + C r||(w,F)||i.||Ff /,^,,,dr; 
^ Jo Jo 

<^*ll(Vw,VF)||i.dr 



w VF : Fdxdr 
ft 



E'^F : Vwdxdr 



+ C / ||(w,F)|li.||Fr._, 



+ n/p + 2/q 



dr. 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



(5.30) 



where (E> denotes the tensor product. Returning back to the estimate (|5.26l) and putting 
(E^-dOOl) together yield that 



||(w(t),F(t))||i.+2 / ||(Vw(r),VF(r))||i.dr< ||(wo,Fo)||i. 



+c f\\Mh + \\E\\h)i\\ur 

Jo ^p. 



l + n/p + 2/g 



+ 11^111-. 



+ n/p + 2/q 



)dr. 



(5.31) 
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This estimate together with Gronwall's inequahty yield the desired estimate (|5.3p immedi- 
ately. We complete the proof of Proposition [5?2l ^ 
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